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Coherent Diffraction Imaging of Single 95nm Nanowires 
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Photonic or electronic confinement effects in nanostructures become significant when one of their 
dimension is in the 5-300 nm range. Improving their development requires the ability to study their 
structure - shape, strain field, interdiffusion maps - using novel techniques. We have used coherent 
diffraction imaging to record the 3-dimensionnal scattered intensity of single silicon nanowires with 
a lateral size smaller than 100 nm. We show that this intensity can be used to recover the hexagonal 
shape of the nanowire with a 28 nm resolution. The article also discusses limits of the method in 
terms of radiation damage. 
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Vertical semiconductor nanowires (NWs) are devel- 
oped as new, nanoscale building blocks for future elec- 
tronic and photonic devices. For all applications, the 
properties depend strongly on the individual crystal 
structure and the average characteristics of the assem- 
bly of NWs. For some applications (resonant sensors, 
optical microcavities), device operation requires precise 
and matching characteristics (like diameters, longitudi- 
nal insertions or core/shell thicknesses) on many NWs. 
To characterize the structural properties of these new 
materials, transmission electron microscopy (TEM) is a 
standard tool that gives accurate information about indi- 
vidual objects when they are small (diameter < 20 nm). 

X-ray diffraction is also a choice method to study 
nanostructures ji particularly when chemical sensitivity 
(using anomalous scattering 2 ) or strain mapping is re- 
quired. But there is a limited number of X-ray diffrac- 
tion studies on nanowires 3,4,5 , due to the fact that for 
many samples it is not possible to perform a quantita- 
tive study on an assembly of NWs: lack of an epitaxial 
relationship between the substrate and the vertical NWs, 
too large angular distribution of the wires, etc.. In this 
case it is required to conduct diffraction experiments on 
single NWs. To this end, Coherent Diffraction Imaging 
(CDI) has been developed during the last few years £ in 
this method, scattering from a single object is recorded 
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FIG. 1: Scanning electron microscopy images of a sili- 
con nanowire, with an hexagonal-shaped section. The gold 
droplet used as a catalyst for the synthesis can be seen at the 
top of the wire (upper view). The bottom view shows the 
facets and nano- facetting for the same wire. 



around one Bragg peak, and the recorded 3D data can 
be inverted to recover the shape of the scattering object. 
State of the art recent results show that it is now possible 
to study objects down to a few hundreds of nm in size, 
and more importantly that CDI allows to probe defor- 
mation inside a nanocrystal^ Existing studies have been 
focalised mostly on model materials with heavy atoms to 
yield a strong scattering of the X-ray beam (e.g. Pb and 
Au droplets). 

In the case of NWs, optical confinement effects be- 
gin to occur around 100-300 nm diameter (depending on 
the material), and electronic confinement around 5-30 
nm^ It is therefore of particular importance to study 
small wires, with a diameter of the order of 100 nm. To 
study the possibility of using CDI for the study of sub- 
100 nm wires we conducted an experiment with homoge- 
neous nanowires on the ID01 beamline of the European 
Synchrotron Research facility. 

The silicon NWs (see Fig. [1]) were grown by low pres- 
sure (20 mBar) chemical vapor deposition at 650°C via 
the gold-assisted Vapor-Liquid- Solid mechanism on (111) 
silicon substrate. Au droplets defining the NW size were 
obtained by dewetting of 2D layers and silane (15 seem) 
diluted in hydrogen was used as the reactive gas. The ad- 
vantage of these Si NWs is that they are known to grow 
with little or no defects -particularly stacking faults along 
the [111] direction- which is vital in a CDI experiment 
since it is required that all atoms diffract in a coherent 
manner. In the case of faulted NWs the scattered signal 
would be vastly different £ 

During the experiment, the epitaxial wires were re- 
moved from the substrate and deposited on another sil- 
icon substrate, thus producing a random orientation for 
all the wires - this method decreases the likelihood that 
two different wires will diffract simultaneously on the de- 
tector for a given orientation of the sample, so that the 
diffraction of a single wire can be recorded. The sample 
was placed under He atmosphere to prevent oxidizing. 

The experiment was carried out using an X-ray pho- 
ton energy equal to 10 keV, using beryllium compound 
refractive lenses 9 to focus the beam to 8.2 x 13.6 \im 2 
(larger than the NW lengths), with a photon flux equal to 




FIG. 2: (color online) 3D coherent diffraction image of a single 
Si < 111 > nanowire, shown as a projection of multiple semi- 
transparent iso-surface layers, with the intensity increasing 
logarithmically from blue to green and red. The hexagonal 
symmetry of the wire is visible in the 3D diffraction image. 
The scale is given for k = 2 sin 0/ A. 

4x 10 10 ph/s. The scattered intensity was measured using 
a Roper Scientific (SX-1300B) direct-illumination CCD, 
to obtain a maximum resolution and photon-counting 
emciencyJ^ This detector presents 20 fim pixels and was 
placed at 743 mm from the sample. 

In order to collect the 3D scattered intensity around 
one Si (111) Bragg peak, the detector was placed at the 
Bragg angle corresponding to the reflection, and the sam- 
ple was then translated and/or rotated until one NW 
diffracted on the detector. Goniometer translations were 
used to ensure that the NW remains at the intersection of 
the X-ray beam and of the horizontal rotation axis. The 
complete scattering was then recorded by rotating the 
sample over a 1° angular range with 0.01° steps to yield 
a 3D pattern, which is then projected to an orthonormal 
frame of reference in reciprocal space (see Ref. 11 for de- 
tails on the conversion between detector and reciprocal 
space coordinates) ; a view of the corrected 3D pattern 
is shown in Fig. [2j 

In the case of a single, non-strained nano-object sub- 
jected to an incoming coherent plane wave, the scattered 
amplitude A is equal to : 

A(k) = F(k)fT[fi(r)] (1) 

where k = kf — lq is the scattering vector near the 
kBragg Bragg position, F(k) is the structure factor of 
a single unit cell, FT[Q(r)] is the Fourier transform 
(FT) of the shape function fi(r) of the nano-crystal, 
i.e. a function which is equal to 1 inside the crystal 
and outside. 20 As the structure factor is slowly vary- 
ing near an existing Bragg peak, we can further assume 
that F(k) « ^ (keragg) and therefore that the scattered 
amplitude is the 3D FT of the nanocrystal's shape - or 
equivalent ly to its electronic density profile. 

During a CDI experiment only the square modulus of 
the amplitude is collected, and all phase information is 
lost. It is however possible to recover the lost phases 
using a priori information on the object - in our case 
the inverse FT of the scattered amplitude corresponds 
to the projection of the electronic density of the sam- 
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FIG. 3: (color online) Simulation of the truncation effect on 
the reconstruction of a NW shape, (a) original hexagonal 
shape of the NW ; (b) FT of the NW shape: as the original 
shape is cent rosy mmetric, the FT is real, here depicted with 
positive regions in red and negative in blue ; (c) truncation of 
the FT and (d) resulting NW shape calculated by inverse FT, 
in which dips (representing & 15% of the average density) can 
be clearly seen due to the truncation. 



pie and must be real, positive, compact and finite- sized. 
These criterion have been used in a number of iterative 
phase retrieval algorithms— d&lldiL to recover the phase 
corresponding to each point of the 3D reciprocal space. 

In our case the recorded intensity, shown in Fig. [2j ex- 
hibits a clear sixfold symmetry which suggests that the 
original < 111 > wire was hexagonaly-shaped, as ex- 
hibited on most wires by scanning electron microscopy 
(SEM). Furthermore, the calculation of the FT of an 
hexagon shape reveals that the phase is either or tt 
(due to the presence of a center of symmetry), as shown 
in Fig. 02 and presents concentric hexagonal-shaped rings 
with alternated 0/tt phases. In order to reconstruct the 
shape of the NW 2D cross-section^, we projected the 
recorded intensity onto the plane perpendicular to the 
wire and applied an alternated sign (+1/ — 1) to the con- 
centric rings (see Fig. |4k) and then computed the inverse 
FT to yield the shape of the NW, which can be seen in 
Fig. Sb. 

The obtained shape of the NW exhibits the expected 
hexagonal symmetry. In order to check that the 'alter- 
nating 0/tt phases' was the correct solution to the lost 
phases, we checked (see Fig. Hh) that the real part of the 
inverse FT was strictly positive inside the object, and 
much larger than the imaginary part (ratio > 50), 

which validates the phase choice. Note that this outcome 
was not guaranteed, as the original data was not strictly 
centrosy mmetric, and no attempt (e.g. by symmetry av- 
eraging) was made to modify this. 

The size (FWHM distance between opposite faces) of 
the wire is « 95 nm, with a resolution of about 28 nm 
(calculated from Ak max — 0.035nm _1 ) for the recon- 
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FIG. 4: (color online) The real shape (the projection of the electronic density) of the NW can be reconstructed from the 
measured scattered amplitude by recovering the lost phase information, either by predicting the phase values (a-b-c) or by 
using iterative phase retrieval algorithms (d-e-f) (see text for details) (a) projection of the experimental scattered amplitude, 
with model phases - the saturation corresponds to the logarithm of the intensity, and the color to the phase. Relative coordinates 
around the (111) reflection are given in k = 2 sin 0/ A units, (b) NW density cross-section, obtained by inverse Fourier transform 
of (a) ; the distance between opposite facets is equal to w 95 nm. (c) density profiles along the horizontal (p x ) and vertical 
(p y ) directions, which demonstrate (the real components being positive and the imaginary one negligeable) that the choice of 
alternating (+) and (-) sign for successive rings is correct. The dips of the density inside the NW are due to the limited extent 
of the intensity scattered by the silicon NW.(d) Experimental amplitude with the phase recovered using ab initio Fourier- 
recycling algorithms, (e) the best reconstruction obtained and (f) the corresponding cross-sections (positivity is imposed in the 
algorithm) . 



struction ; the NW appears to be slightly elongated in 
one direction, which may be due to an asymmetric NW 
growth, probably due to the original seed combined with 
kinetic limitations. In the extracted density map shown 
in Fig. [4)3, we observe a variation amounting up to 20% 
from the average value. This is very likely due to the 
limited extent of the recorded intensity: only three in- 
terference orders (0,1,2) have been recorded, due to the 
small size of the NW and the weak scattering power of 
the silicon atoms, which both limits the real-space res- 
olution and can create ripples inside the reconstructed 
density ; a simulation of the effect of truncated data is 
presented in Fig. [3] and has allowed us to reproduce the 
electronic density dips (Fig. (3ji). Another source for the 
limited resolution on the sides of the NW comes from the 
zig-zag nature of the NW facets (see Fig. [I]), which blurs 
the average density on the NW borders. 

In order to check that the reconstructed structure us- 
ing alternating phases was the best possible, we also per- 
formed an ab initio reconstruction using a combination 
of the error-reduction, hybrid input/output and charge 
flipping algorithms, using finite support and positivity 
as constraints, following the process depicted in Ref^. 
The best set of phases that we obtained 22 (see Fig.0]i,e,f) 



has an R-factor 12 equal to 0.7%, features the expected al- 
ternating phases, and the reconstructed NW shape also 
presents a quasi-hexagonal shape, thus confirming the 
results of the inversion using "predicted phases". 

One expected problem when studying single, nano- 
sized objects using an X-ray beam is radiation dam- 
age: this is a well known issue for macromolecular 
compound ^ 16 ! 17 , with a dose limit estimated to 2. 10 7 Gy 
(IGy = U/kg). In the case of inorganic compounds 
there is no acknowledged dose limit, as the compounds 
do not present the specific weak bonds (— CO2, S — S) 
that are likely to break under irradiation. 

We have however observed that the NW can break un- 
der the X-ray beam, as is illustrated in Fig. [5l We mea- 
sured the scattering from a single wire by rotating the 
sample over a 0.8° angular range with 0.02° steps, and 
a 50 s exposure time per image. This was repeated four 
times in order to accumulate more statistics, and the evo- 
lution of the projection of the 3D recorded scattering is 
shown in Fig. [5j During this measurement, the flux of 
the experiment was equal to « 4ph/s/A 2 . Given the sil- 
icon absorption cross-section (JE=iokeV — 1-5 x 10 _5 A 2 , 
this results in a absorbed power per atom of 0.6 eV/s, or 
equivalently 2 x 10 6 Gy/s, i.e. a total dose of 8 x 10 9 Gy 
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FIG. 5: (color online) Effect of radiation damage: (a,b,c) evo- 
lution of the 3D diffraction pattern of a single Si NW around 
the (111) reflection for successive scans, with a total exposure 
time of 4000 s per scan (a fourth scan is not presented here). 
The projection is presented in the direction perpendicular to 
the NW axis, using a logarithmic colour scheme. The short 
vertical oscillations are due to the finite length of the wire, 
and almost disappear in later scans, due to the damage to the 
NW. The fringes which can be seen parallel to the diagonal 
of the image correspond to the FT of the NW cross- sect ion, 
as seen in Fig. [2] for another wire. The split of the main 
diffraction peak could be due to the breaking of the NW: (d) 
simulated scattered amplitude around the (111) reflection of 
a broken wire, with two parts with an equal length of 400 nm 
along < 111 > and a width of 105 nm along < 110 >, with a 
separation between the two parts equal to half the interpla- 
nar distance dm. This separation leads to two maxima at 
the (111) reflection position. 



per scan. 

During the experiment, the absorbed energy is not ac- 
cumulated and can be evacuated through fluorescence, 
thermal radiation or conductivity, but point defects will 



nevertheless be created in the sample. The presence of 
already existing defects (such as a stacking fault -a com- 
mon occurence in NW— , point defects or dislocations) 
can however be considered as a "weak" part of the struc- 
ture, which could lead to a broken wire as exhibited in 
Fig. [5l However new experiments (e.g. with a more fo- 
cused and intense beam) would be necessary to quantify 
the relationship between the dose, existing faults, and 
the effects of radiation damage to a NW. 

In conclusion, we have demonstrated that it was pos- 
sible to use coherent diffraction imaging on single ~ 
100 nm NW to recover their shape with a ~ 28 nm res- 
olution, even in the case of weak scatterers like silicon. 
Evidently the interest of this result does not lie in the 
reconstruction of the shape of homogeneous NW, since 
scanning and transmission electron microscopy routinely 
yield the external shape of nano-sized object with a bet- 
ter resolution. However this technique will be partic- 
ularly useful in the case of heterogeneous NW, i.e. 
with longitudinal or radial (core-shell) heterostructures: 
these type of NW exhibit strain fields and chemical gradi- 
ents at the interface for which X-ray diffraction remains 
the most sensitive technique, and will be the focus of up- 
coming experiments. Such an experiment would allow to 
recover the 3D complex field describing both the density 
and the distortion of the lattice, as depicted in Ref 6 . 

Moreover constant progress is made for X-ray optics in 
order to improve the coherence and flux on synchrotron 
beamlines - e.g. as of the writing of this article the fo- 
cused beam size on the ESRF ID01 beamline is much 
smaller (10 10 ph/s in 2.5 x 0.5 jam 2 ), yielding a flux 20 
times more intense than during our experiment, which 
should allow to study smaller NWs with a better resolu- 
tion. 
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We have also studied the case of GaAs/GaP [111] NWs 
which are known for their stacking faults along the growth 
direction. In such a case, the CDI signal recorded is essen- 



tially be a function of the faults nature and distribution 
along the NW axis, instead of the shape and elastic defor- 
mation of the wire. 

Note that this approach is only correct as long as the crys- 
tal consists of complete unit cells - a reasonable approx- 
imation as long as the dimensions of the wire are much 
larger than the unit cell parameters. 

As the NW is much longer (several \im vs. « 100 nra), the 
FT of its shape is essentially 2D. 

The optimization used 200 to 400 cycles, with a fixed 
square support more than twice the size of the NW, and 
was repeated 400 times to avoid stagnation. 



